function quadRuleTr = quadRule2D()

quadRuleTr{3}.nodes =  [0.659027622374092    0.10903900907287759   0.23193336855303046 ...
  0.23193336855303046   0.10903900907287759  0.659027622374092;  ...
  0.23193336855303046   0.659027622374092    0.10903900907287759 ...
  0.659027622374092     0.23193336855303046  0.10903900907287759;...
  0.10903900907287759   0.23193336855303046  0.659027622374092     ...
  0.10903900907287759   0.659027622374092    0.23193336855303046];

quadRuleTr{3}.weights = ...
  [0.16666666666666666 0.16666666666666666 0.16666666666666666 ...
   0.16666666666666666 0.16666666666666666 0.16666666666666666];

quadRuleTr{5}.nodes =  [0.3333333333333333  0.10128650732345634   0.4701420641051151   ...
  0.10128650732345634   0.4701420641051151  0.7974269853530873    0.05971587178976982; ...
  0.3333333333333333    0.10128650732345634 0.4701420641051151  ...
  0.7974269853530873    0.05971587178976982 0.10128650732345634   0.4701420641051151;  ...
  0.3333333333333333    0.7974269853530873  0.05971587178976982 ...
  0.10128650732345634   0.4701420641051151  0.10128650732345634   0.4701420641051151];

quadRuleTr{5}.weights = ...
  [0.225                 0.12593918054482717   0.1323941527885062  ...
   0.12593918054482717   0.1323941527885062    0.12593918054482717 ...
   0.1323941527885062 ];

quadRuleTr{7}.nodes = [0.06493051315916486   0.06493051315916486   0.8701389736816703 ...
  0.2845755842491703   0.3135591843849315    0.19838447668150672 ...
  0.642577343822696    0.517039939069323     0.04386347179237249 ...
  0.517039939069323    0.04386347179237249   0.19838447668150672 ...
  0.642577343822696    0.2845755842491703    0.3135591843849315; ...
  0.06493051315916486  0.8701389736816703    0.06493051315916486 ...
  0.517039939069323    0.04386347179237249   0.2845755842491703  ...
  0.3135591843849315   0.19838447668150672   0.642577343822696   ...
  0.2845755842491703   0.3135591843849315    0.517039939069323   ...
  0.04386347179237249  0.19838447668150672   0.642577343822696;  ...
  0.8701389736816703   0.06493051315916486   0.06493051315916486 ...
  0.19838447668150672  0.642577343822696     0.517039939069323   ...
  0.04386347179237249  0.2845755842491703    0.3135591843849315  ...
  0.19838447668150672  0.642577343822696     0.2845755842491703  ...
  0.3135591843849315   0.517039939069323     0.04386347179237249];

quadRuleTr{7}.weights = [0.053077801790232484 0.053077801790232484 0.053077801790232484 ...
   0.07085308369213365   0.069274682079417    0.07085308369213365  ...
   0.069274682079417     0.07085308369213365  0.069274682079417    ...
   0.07085308369213365   0.069274682079417    0.07085308369213365  ...
   0.069274682079417     0.07085308369213365  0.069274682079417];
end